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SUMMARY 

Equations are reviewed which show that at equilibrium fugaclty and con- 
centration gradients can exist In gravitational fields. At equilibrium, the 
logarithm of the ratio of the fugacltles of a species at two different loca- 
tions In a gravitational field Is proportional to the difference In the heights 
of the two locations and the molecular weight of the species. An analogous 
relation holds for the concentration ratios In a multicomponent system. The 
ratio Is calculated for a variety of examples. The kinetics for the general 
process are derived, and the time required to approach equilibrium Is calcu- 
lated for several systems. The following special topics are discussed: Ionic 

solutions, polymers, multiphase systems, hydrostatic pressure, osmotic pres- 
sure, and solubility gradients In a gravity field. 


INTRODUCTION 

Application of the relationship that exists between fugaclty and position 
In a gravity field might lead to a better understanding of the behavior of 
both living and nonliving matter In a gravity field. For Instance, this rela- 
tionship might be Involved In the mechanism by which the roots of a plant are 
directed downward. Eventually, the resultant knowledge may suggest micro- 
gravity experiments that will add still more to the understanding of gravity. 

Consider a material In a gravitational field. The upper portion of the 
material has a higher gravitational potential energy than the lower portion. 

Of course, the difference Is small compared with most other effects. In zero 
gravity this effect does not exist, and all parts of a material have the same 
potential energy. The existence of a gravitational potential gradient can 
affect the behavior of materials. 

Sometimes, this fact Is not recognized even though the effect was treated 
years ago by Guggenheim (ref. 1) and others. Guggenheim began with the fol- 
lowing equation. 


d(v^ + M^ q>) = 0 


( 1 ) 


where p Is the chemical potential; <p Is the gravitational potential; and 
1 refers to the 1'th species In a solution. Equation (1) Ignores contribu- 
tions from surfaces and fields other than gravitational. Via equation (1) 
Guggenheim expressed the change In mole fraction, N^, of a species In an 
Ideal solution as a function of the change In gravitational potential. 
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where M refers to molecular weight; V, to partial molal volume; X(NkMfc) Is 
the effective molecular weight of a mole of the solution; and X(N|<V|<) Is the 
volume of a mole of solution. The ratio, X(N| c Mk)/X(N| < Vk) , equals the density 
of solution, p. Koenig (ref. 2) extended Guggenheim's equation to apply to 
nonideal solution. More recently Guggenheim (ref. 3) has expressed the fugac- 
1 ty of species In a mixture of perfect gases as a function of position In a 
gravitational field. 
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where p Is the fugaclty (Idealized pressure); and M, the molecular weight. 
Guggenheim's equations are difficult to apply In the form they are given. 
Because a single phase Is sometimes defined as a system of uniform composi- 
tion, and a gravitation field requires a nonuniform composition at equilib- 
rium, Guggenheim refers to an equilibrium system In a gravitational field as a 
multiphase system. Thus, by his definition, an equilibrium aqueous solution 
of sodium chloride In the earth gravitational field should be referred to as a 
multiphase system. Such terminology Is confusing and, therefore. Is not used 
In this report. 


In the case of an Ideal gas, the equations that have been derived are In 
a more usable form (ref. 4) 


In 



-Mg(h - h Q ) 
RT 


(4) 


where P Is the pressure of the gas at the higher altitude h; P 0 Is the 
pressure at the altitude h 0 ; M Is the molecular weight; g Is the accelera- 
tion of gravity; R Is the gas constant; and T Is the absolute temperature. 
Equation (4) can be derived In several ways. Independently of Guggenheim's 
equations. In a mixture of gases, each component obeys equation (4) Individu- 
ally where M, P, and P 0 refer to the same component. The concentration 
gradient In the earth's atmosphere Is due at least In part to the relationship 
expressed In equation (4). Equation (4) also applies to the separation of 
Isotopic gaseous species under the high gravity conditions of a centrifuge. 


In 1908, Perrin (ref. 5) modified equation (4) In order to apply It to a 
suspension of particles In a liquid. 


IP m ^ | nr 3 gh(p - p') (5) 

where r Is the radius of the particles; N Is Avogadro's number; p l Is the 
density of a particle; p Is that of the liquid; n and n 0 refer to the 
concentration of the particles at two different heights; and h Is the dif- 
ference In their heights. From experimental data and equation (5), Perrin 
calculated Avogadro's number to about a 15 percent accuracy. His experiments 
validated equation (5). Since then, however. It has been shown that 
equation (5) holds only In a relative narrow range of h (ref. 6). 
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Block et al. (ref. 7) have found concentration gradients to occur In 
gravitational fields In excess of 10$ times normal gravity for a liquid 
solution several degrees above Its liquid - liquid critical solution point. 
The solution was a mixture of n-hexane and perf luoro-n-hexane. Similarly, 
Moldover et al. (ref. 8) have reported density gradients to exist In gravity 
fields for systems near their gas - liquid critical point. Other Investi- 
gators have shown similar results. 

Diffusion kinetics for the attainment of concentration equilibrium In a 
gravitational field appears to be lacking from the literature although diffu- 
sion of Ions In an electrical field has been treated (ref. 9). 

The purpose of this report Is to present the appropriate equilibrium 
equations In usable forms, to apply them first to simple examples, to derive 
the applicable kinetic equations, to evaluate the time required for approach 
to equilibrium, and finally to apply the equilibrium equations to certain 
special cases such as Ionic solutions and osmotic pressure. 


EQUILIBRIUM EQUATIONS AND CALCULATIONS 

In the first part of this report, our attention Is directed toward only 
single phases In a uniform gravitational field. We define h-| as a distance 
between the gravitational source and some position In the phase and h 2 as 
the distance between another position In phase and the gravitational source. 

If the difference between these values, h, Is small with respect to the values 
themselves, the gravitational potential may be expressed as a linear function 
of h. In differential form the relationship becomes: 


d<p = - g dh (6) 

where g Is the acceleration of gravity, and h Is measured toward the grav- 
itational source. 


Concentration 

To make equation (2) more usable In predicting concentration gradients In 
the condensed state, equation (6) Is used to express the gravitational poten- 
tial difference. Combining equations (2) and (6) yields: 

RT dN, 

RTd £n(N i ) = — j^- 1 = (M 1 - V^) g dh (7) 

If J refers to the solvent, dNj = - I dN^ . Combining this equation 
with equation (7), 

RT dNj = - |lN 1 (M 1 - V lP )}g dh (8) 

Moreover, the average molecular weight of the solution, M = NjMj + X(N^M^); 
the volume of a mole of solution, V = N j V j + X(N^); and M = V P . Thus 
X(N^(M^ - V 1P )} = -Nj(Mj - Vjp) . Combining this equation with equation (8), 
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RT dN. 

-jj-l . (Mj - V jP )g dh (9) 

Equation (9) states that If all the solutes obey equation (7), the solvent 
obeys the same relationship. 

Equation (7) can be Integrated If p Is assumed to be constant. This 
assumption Is only valid for a dilute noncompresslble solution. 

/N"\ 

RT * n UjJ = g h ( H i - V> - gh(V 1 /> 1 - V iP ) = ghB 1 (10) 

where Is the partial molal volume of species 1 In solution; Is the 

molar volume of pure species 1; and p^ Is Its density. Equation (10) Is 
completely valid for the solute In a dilute Ideal noncompresslble solutions 
when the distance to the gravity source Is large. The product, V^p, Is the 
weight of solution corresponding to a molal volume of species 1. The mass 
difference, p, corresponds to a buoyancy factor, . When a species 

has a positive buoyancy factor, the species can be said to be more dense than 
the solution so that at equilibrium the concentration (mole fraction) of the 
species will be greater at the bottom than at the top. When a solute has a 
negative buoyancy factor. Its concentration at equilibrium will be greater at 
the top than at the bottom. 

The term Ideal system (or perfectly Ideal solution) Is generally taken to 
mean that both the solvent and solute behave Ideally; Raoult's law Is obeyed. 

In this case, Is equal to the molar volume of the pure species because 
V 1 will not be a function of concentration. Thus, 

RT ,n (s|)‘ H i( 1 ' ^) 9h (11) 

where 1 can refer again only to a solute. The examples of truly Ideal sys- 
tems are rare. 

In general, one must be satisfied with the treatment of "Ideal solutions" 
In which the solutes obey only Henry's law. Henry's law requires that 
does not vary with concentration In the low concentration region. Outside 
this region equation (10) should be considered to be only an approximation; 
the value to be used for outside the low concentration region should 
likely be the value of for the concentration of Interest. 

It should be noted that equation (10) does not predict a constant concen- 
tration gradient throughout a column of material In a gravitational field. As 
a result the equilibrium concentration at midheight. In general. Is not equal 
to the Initial concentration, Ni f0 . The value of the midheight concentration, 
Ni m . Is dependent upon the variation of the cross section of the column with 
height as well as on the other gravitational parameters. To derive an expres- 
sion for N* m , we have assumed a uniform cross section. In a dilute solu- 
tion, the average concentration which Is equal to N^ >0 may be set equal to 
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X(N^ dX)/h, where X Is the distance measured from the top, and Is the 
corresponding concentration. The variable, N^, may be replaced by exp{(M^ 

- V^p)gX/RT}, where Is the mole fraction at the top. Integration, 
rearrangement, and substitution yield: 


_Nj_ (H 1 - V iP )gh 

N 1 o " T (M 1 - ^P)flh 
’ RT exp - J RT 1 - 1 



( 12 ) 


From equation (10) It can be shown that m equals sqr(N^' 
duclng this Into equation (12) yields: 


Nj). Intro- 


N 

N 


1 ,m 
1,0 



(13) 


In figure 1, the ratio N^ >m /N^ >0 Is plotted as a function of The 

midheight concentration Is’less’than the Initial concentration, both when the 
ratio N^'/n] Is greater than unity (positive buoyancy factor) and when It Is 
less than unity (negative buoyancy factor). When N^'/Nj approaches unity 
(corresponding to a small product of g, h, and the buoyancy factor), N^ >m /N^ >0 
also approaches unity. The fact that m Is not equal to N^ 0 for a solu- 
tion In a gravitational field means that Hs center of mass Is not at Its geo- 
metric center. 


Application to Salts 

A question arises when applying the equations just derived to Ionized 
salts. Should one treat the salt or the Ion as the solute species? 

Equation (1) which Is the starting point of the derivation Is for a single 
species. However, when the species Is an Ion, It must be modified to Include 
and the electrical potential, Z^: d(vi + p] + Z^) = 0. To derive a usable 

equation from this modified equation, Z^ must be eliminated by adding 
together the equations for all the species In a molecule of the salt. From 
such a composite equation, an equation similar to equation (7) can be derived. 

RT dN (M - V p)g dh 

RTd tn(N^ ) = — ^ (14) 

where Y Is the total number of Ions present In a molecule; s refers to the 
salt; and 1 can refer to either the cation or the anion. The 1 can also 
refer to the salt In the sense that the concentration of the 'Ionized' salt 
must vary with the concentration of the anion and the cation. Remember, how- 
ever, that the concentration of the 'nonlonlzed' salt species Itself Is zero 
because we assumed a completely Ionized salt. However If some of the salt 
were nonlonlzed, that component of the salt would obey equation (10). 
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Equation (14) applies only when a single salt Is present In the solution. The 
problem of the presence of more than one salt Is treated In a later section. 

By the comparison of equations (10) and (14) one sees that the buoyancy factor 
for an Ionized salt, B s , Is (M s - 9 s p)/Y. Thus equation (10) Is applicable to 
an Ionized salt If the proper buoyancy factor Is used. 


Fugaclty 

A gravitation field also gives rise to a fugaclty gradient In a column of 
material. This gradient may be looked upon as a combination of two effects: 
the concentration gradient just discussed and the hydrostatic pressure of the 
column Itself. The effect of pressure on fugaclty Is well known (ref. 10). 


RT 



or more precisely. 


RT 


d 

dP 



(15) 


where f^ Is the fugaclty of species 1 without the applied pressure, P; 

f}' Is Its fugaclty with the applied pressure; and Is the partial 
molal volume of the species. Because the hydrostatic pressure Is equal pgh, 
equation (15) can be written as follows: 


RT 


d * n ( f l. 0 > 

dh 


v^g 


( 16 ) 


where the subscript o Indicates the condition of homogeneous composition. 
Thus even without a concentration gradient, a fugaclty gradient exists In a 
column of material In a gravitational field. This equation applies to both 
the major and minor components (solvents and solutes). 


Because Henry's law Is assumed, the ratio, f^/N^, Is a constant for 
an nonlonlzed material at any specified position In the column, and f ^ /N^ 

= f 1 0 /Ni o» fl being the equilibrium value of the fugaclty. For an Ionized 
salt’ f s /l^i Y = f s ,o/Nl,o Y * Because N^ >0 Is not a function of h, one 
may write 

d an(f $ ) d Jtn( f s 0 ) Yd %n(HJ (17) 

dh = dlT* + dh 


For an nonlonlzed material, the s subscripts become 1 ' s and the Y becomes 
one. Combining equations (7), (16), and (17) yields 

d JLn( f . ) 

RT - H lfl 
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or 


RT fcn^jj = M^h (18) 

where is the equilibrium fugaclty value for species 1 at the top of 
the column, and f\ Is Its equilibrium value at the bottom. The Integration 
step requires g to be constant. Equation (18) holds for liquids, solids 
(Including Ionized salts), gases, solutes and even solvents. Moreover, because 
equation (18) does not specify properties of the solution. It holds even though 
Henry's law Is not obeyed. 

Equation (18) as applied to gases may also be derived by combining equa- 
tions (3) and (6). For an Ideal gas, one may also substitute pressures for 
f ugacltles and thus derive equation (4). Thus, In a mixture of Ideal gases, 
equation (4) will apply Individually to each component. Of course, when a gas 
or gas mixture deviates from Ideality, equation (4) can still be used as an 
approximation. 


Numerical Evaluation 

Calculation of a concentration gradient expected In a gravitational field 
at equilibrium requires the knowledge of of the solute. This value Is 

calculable from the variation of density of a solution with concentration. 


(100 - pj) (100 - P p 



where p' Is the density of a solution containing <p\ u weight percent of 
solute 1; p", the density of a solution containing weight percent of 
the solute; and M^, the molecular weight of the solute. Equation (19) Is 
approximate but becomes precise In the limit where p^ approaches p\. 

The required density data as a function of concentration are available 
from a variety of sources (e.g., aqueous solutions from ref. 11). Via such 
data and equation (19), we have calculated the values of and thereby 

the buoyancy factors for a variety of solutes In water. These are listed In 
table I. The molecular weight used for a solute Is that of Its anhydrous form. 
The value of Y (number of species per molecule) that was used Is listed In the 
second column. For nonionic materials, Y Is one. For sodium sulfate which Is 
a completely Ionized salt, Y Is three (two sodium Ions and one sulfate Ion). 
The concentration ratios given for It In table I are therefore those for the 
Ions. For sodium chloride, also completely Ionized, Y Is two (one sodium Ion 
and one chloride Ion). For only slightly Ionized materials such as acetic 
acid, ammonium hydroxide, and phosphoric acid, the value of Y given In 
table I Is one; therefore, the concentration ratios given refer to those of 
the nonlonlzed species which Is the major component. Trichloroacetic acid Is 
somewhat of a special case because It Is almost but not completely Ionized; In 
this case, the value used for Y Is two. It Is Interesting to note that at 
Infinite dilution, a partially Ionized material would become completely 
Ionized, and the Y to be used would be that of the completely Ionized 
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material. Under such consideration, the value of Y for acetic acid would be 
two, and the concentration ratios would refer to those of the Ions. However, 
the Y's used In table I are ones of practical Interest, those for 'reasonably 
dilute solutions for which such materials are essentially not Ionized. 


Selection of Parameters 

At equilibrium, concentration gradients are not found under normal gravity 
conditions. Only when large values of h (ref. 4) or large values of g are 
Involved (refs. 7 and 8), have concentration gradients been observed. The 
three sets of parameters used In this report were selected In order to encom- 
pass the conditions where the gravitational effect Is expected to be detected 
and where It Is not. Parameters 'a' are 1 G (l.e., g = 980 cm/secO and a 
1 m high column; parameters 1 b ‘ are 100 000 G's and a 1 cm high column; and 
parameters * c 1 are a 10 000 m column at 1 G. This height approximates the 
depth of the 5 mile ocean. Because the gravitational concentration effect Is 
dependent on the product of g and h (eq. (10)), ratio c Is also valid for 
one million G's and h = 1 cm. In table I the ratios of the lower to upper 
concentrations are given for these three sets of gravitational parameters. 
Equation (10) was used to make the calculations. 

Examination of table I shows that the concentration ratios 'a' (corre- 
sponding to parameters a) are generally very close to unity; this Is In 
accord with observations under normal conditions (a relatively small value for 
the product of g and h.) Calculated concentration differences between the 
top and bottom are In most cases less than one part per thousand. Only In the 
case of high molecular weight material (e.g., dextran) does the calculated 
ratio (1.122) depart significantly from unity. 


Ultracentrifuge Conditions 

In contrast to concentration ratios close to unity, table I shows that a 
considerable number of concentration ratios 'b' (for 100 000 G's and a 1 cm 
high column of solution: ultracentrifuge conditions) are greater than 1.1. 

This corresponds to concentration differences from top to bottom of more than 
10 percent. Thus, concentration gradients of many dissolved materials should 
be measurable under ultracentrifuge conditions. There are several extreme 
examples In table I. One Is the ratio for lead nitrate, 307, which means a 
concentration gradient of over 30 000 percent. It would appear from table I 
that the ratio Is greatest for high molecular weight materials, e.g., poly- 
mers. The ratio for dextran Is greater than 10 38 ; the ratio for Inulln, 

2767 As discussed earlier, when the buoyancy factor Is negative, the concen- 
tration of the solute at the top of a column Is greater than at the bottom. 
This Is demonstrated In the case of the organic liquids: acetone, ethanol, 

and methanol. Ammonium hydroxide also has a negative buoyancy factor. 


Effect In Ocean 

The concentration ratios * c * In table I (for 1 G and a 10 000 m deep 
ocean) deviate even more from unity than the concentration ratios b. Many of 
the concentration ratios ' c ' are greater than 10. This would Imply that at 
equilibrium, many dissolved materials should tend to concentrate at the bottom 
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of the ocean and lead to the purification of Its surface regions. For 

Instance, sugars such as glucose would be about 15 times more concentrated at 

the ocean's bottom than at Its surface. The ratio for copper sulfate Is over 
26, meaning that at equilibrium a copper salt dissolved In the ocean water 
would tend to concentrate at the bottom, leaving a much lower concentration 
near the surface. In contrast, the ratio for sodium chloride Is only about 
2.2 which means that the concentration of this salt near the surface would be 

comparable with Its concentration near the bottom. When more than one salt 

exists In solution, equation (10) Is not completely valid because the Ions 
which compose the Individual salts can Interact. This topic Is discussed 
later In more detail . 


Higher Concentrations and Fugacltles 

Also Included In table I Is a 'range' column which Is the percentage range 
In which the partial molal volume of the solute varies only 1 percent. For 
many salts this Is only up to 0.5 or 1.5 percent. However, for many organic 
materials, the range extends to 10 percent. Acetic acid extends to 30 percent; 
fructose and methanol extend to over 50 percent. For this range of concentra- 
tions, the concentration ratios d were calculated assuming a constant average 
density and the parameters of 100 000 G's and h = 1 cm. In general, the ratio 
d does not vary more than about 1 percent from the corresponding ratio b. 

The greatest variation, about 5 percent. Is for fructose; but this Is for the 
concentration range up to 52 percent. The final column In table I gives the 
equilibrium fugaclty ratio for 100 000 G's and a 1 cm high column. 


Nonaqueous Systems 

Equation (10) may, of course, be applied to nonaqueous solutions. This 
has been done In table II using data In reference 12. Table II gives the par- 
tial molal volumes of solutes In the Indicated solvents, the concentration 
ratios for particular values of h and G, and the fugaclty ratios. (The 
fugaclty ratio of a particular species Is Independent of the solvent.) In the 
first part of table II one can find data for organic solvents, while the last 
part gives data for some metals. Again for the condition h = 100 cm and nor- 
mal gravity, the concentration ratios are very close to unity. However, In a 
few cases, the difference between the top and bottom concentrations may exceed 
0.1 percent; for example, the ratio a for tin In gold Is 0.9988. 

In general, ratios b In table II (h = 1 cm and 100 000 G's) just as 
those In table I deviate noticeably from unity. The exceptions are aluminum 
In magnesium and Iron In copper. Again for solutes with high molecular weights 
(e.g., tristearin) the ratios b are unusually large. For many alloys, the 
ratios are less than unity because the corresponding buoyancy factors are 
negative. 

One of the first solutes listed In table II Is sodium nitrate. Its ratio 
data In methylamlne can be compared with Its ratio data In water given In 
table I. Its calculated ratio b Is 1.122 In water while It Is 136 In methyl- 
amine. This large calculated value In methylamlne Is the result of an apparent 
large negative value used for Its partial molal volume. 
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In reference 12, the density of the tristearin - bromocamphor system Is 
given over Its entire composition range. This allowed two sets of calculations 
to be made: one where the tristearin Is assumed to be the solvent and one 

where the bromocamphor Is taken as the solvent. The calculations show bromo- 
camphor to have a positive buoyancy factor when assumed to be the solute, and 
the tristearin to have a negative buoyancy factor when It Is taken as the sol- 
ute. The calculations Imply that regardless of the concentration, the bromo- 
camphor Is expected to be more concentrated on the bottom and tristearin to be 
more concentrated on the top. This brings to mind a metaphysical necessity 
for all two component solutions: If one component Is more concentrated on the 

top, the other must be more concentrated at the bottom. 

The tristearin - bromocamphor system Is useful for another purpose. The 
density data are available for three different temperatures. In table II, 
this data has been used to show that the concentration ratio Is a function of 
temperature as would be expected from the fact that the buoyancy factor varies 
with temperature. 


KINETIC EQUATIONS AND CALCULATIONS 

Although we have seen that at equilibrium, a gravitational field can pro- 
duce a detectable concentration gradients In a variety of systems, one must be 
sure that equilibrium of this sort can be approached In a reasonable time 
frame. 


The equilibrium which we are discussing Is brought about by diffusion. 
However, normal diffusion equations predict a uniform equilibrium concentra- 
tion. It becomes obvious that for our application, the diffusion equations 
must be derived Incorporating the effect of a gravitational field. To treat 
the effect of a gravitational or other type of field, one needs to consider 
only one dimensional diffusion. When a gravitational field Is not considered, 
the following one dimensional equations are applicable (ref. 13). 


F 


1 


dC. 


'1 dX 


( 20 ) 


and 




( 21 ) 


where Is the flow (g/cm 2 /sec); D, the diffusion coefficient (cm 2 /sec); 

C^, the concentration (g/cm 3 ); X, the distance (cm); and t, the time (sec). 
Diffusion proceeds parallel to the concentration gradient. To Include the 
effect of gravity on diffusion, one must return to the derivation of the dif- 
fusion equations. 


Jump Frequency Approach 

One of the derivations for the diffusion equations starts with the concept 
of a molecule jumping from one position to another In a material. Reference 
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to figure 2 will aid In the visualization. One may set up the following equa- 
tion for a solution that obeys Henry's law. 



where Is the flow; «■) , the proportionality constant; X', the distance 

parameter of the starting position (being farther away from the gravitational 
source); X", the distance parameter of the position after the jump (being 
closer to the gravitational source); C^, the concentration of ( species 1 at 
X'; Ci, the concentration at X"; k, the Boltzman constant; , the extra 
energy of the system due to the presence of an I'th molecule, before the jump; 
E? , the corresponding value after the jump; and Q-j , the extra energy of the 
system due to the presence of the I'th molecule In the activated position for 
diffusion (essentially halfway between the Initial and final positions). 


As X" approaches X', equation (22) becomes: 



(23) 


the Ci may be written simply as C^ ; the exp{(-Qi + Ei)/kT} may be 
replaced by Q\, the diffusion coefficient; and i^/RJ may be substituted for 
E^/kT where E^ Is the extra energy of the system due to the presence of a 
mole of species 1. Thus, 



The first term In equation (24) corresponds to the relation expressed In 
equation (20). The second term In equation (24) corresponds to the change In 
potential energy of the system due to a molecule's position In a gravity field 
as the molecule moves toward the gravitational source. Thus far, equation 
(24) Is applicable to any type of field because we have not yet Introduced the 
actual dependence of the energy upon position In a gravitational field. Note 
that equation (24) has the same form as that for an Ion In an electrical field 
(ref. 9, eq. (5) In Chap. 6). 


The change of the concentration with time Is dependent solely on the flow 
gradient as described In equation (21). It Is Independent of any field which 
might exist. Combining equations (21) and (24) yields: 





(25) 


Energy Versus Height 

Equation (26) relates the change of Ei with X within a gravity field. 




( 26 ) 
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where g and are assumed not to be a function of X so that d 2 E^/dX 2 
= 0. In equation (26) as In earlier equations, the expression Is the 

weight of solution corresponding to the volume of one mole of the 1'th species. 
Substituting equation (26) Into equations (24) and (25) yields, respectively: 



Equilibrium 

At equilibrium, In equation (27) becomes zero. The resultant equation 
may be Integrated If the solution Is dilute and noncompresslble. 

RT 9hB 1 (29) 

The c}'/c} may be replaced by N^'/N^, and equation (29) becomes 
equation (10). It should be noted that for the derivation of equation (29), 
the solution did not have to obey Raoult's law; only Henry's law has to be 
obeyed. 


Alternate Derivation 


Equations (27) and (28) can be derived In a simpler and less rigorous 
fashion. In this alternate derivation, one can determine more easily the 
expression to be used for If the diffusing species Is an Ion. The 
derivation Involves subtracting the equilibrium value of dC/dX from Its 
actual value In equation (20). Thus, 



(30) 


where the subscript e refers to the equilibrium condition. For a dilute 
solution, equation (14) can be written as follows: 



(31) 


where the Y = 1 for an uncharged diffusing species. Combining equation (30) 
and (31) result In the diffusion equation for an Ion. 


p D K c a v> | 

F 1 " D 1 |dX " C 1 9 YRT ( 


( 32 ) 
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Equation (32) Is just a more detailed way of writing equation (27). _ Thus, 
equations (27) and (28) apply to Ionic diffusion where B 1 = (M $ - V $ p)/Y. 

Program 

Equation (28) Is not readily Integrated. But with modern computers. 
Integration Is not required to make equation (28) useful. A simple Basic com- 
puter program Implementing equation (28) Is given In table III; It yields the 
concentration and concentration gradient profiles as a function of time. The 
program Is written for an Initially uniform concentration and with the follow- 
ing set of default parameters: buoyancy factor equal 100 g/mol; 100 000 G's; 

D = lxlO -5 cm 2 /sec; h = 1 cm; and 20 finite elements. 


Kinetic Plots 

Our calculations via the program In table III are shown graphically In 
figures 3 to 8. The solid lines In figure 3 gives the concentration profile 
for various periods of time, using the default parameters. These default 
parameters represent a typical solute - liquid solvent system under ultra- 
centrifuge conditions. The points represent calculations using 100 finite 
elements; these are In good agreement with Its respective line for which only 
20 finite elements were used. Note In figure 3 that Initially only the con- 
centrations at the top and bottom of the column change. This Is In accord 
with equation (28): Initially, there are no concentration gradients present 

except for the discontinuities at the top and the bottom. Thus Initially 
according to equation (28), only the concentrations at the top and bottom 
should change. As time proceeds, the concentration gradients move toward the 
middle of the column accompanied by changes In the concentrations. 

One can also see In figure 3, that equilibrium Is approached In about 
50 000 sec (about 13 hr). This Is a reasonable time frame for an ultracentri- 
fuge experiment Intended to measure such a concentration gradient. Figure 4 
Is a similarly calculated profile for the concentration of sodium chloride In 
the ocean, using Its buoyancy factor from table I and a diffusion coefficient 
of 1.54xl0 -5 cm 2 /sec (ref. 14). For this system, equilibrium Is approached In 
about 10 +16 sec (about 500 million years). Because of this long time period 
and the almost constant motion of the ocean, It Is doubtful that gravitation 

equilibrium of sodium chloride In the ocean will ever be approached. 

Figure 5 shows a "concentration" profile that allows a fuller comprehen- 
sion of the approach to equilibrium. In figure 5 the ratio of the concentra- 
tion change to the concentration change at equilibrium Is plotted as a function 
of the distance from the bottom of the column. Thus, the ordinate represents 
the fraction of the equilibrium achieved for the time period Indicated for a 
particular curve. The default parameters were used. A similar type profile 

for the concentration gradient Is given In figure 6. Both figures 5 and 6 

show that equilibrium Is very closely approached after 50 000 sec. Figure 6 
shows that the concentration gradients at the top and bottom are achieved 
Immediately. But this does not mean that the concentration Itself reaches 
equilibrium; see figure 5. 

The Irregularities near the middle of the curves In figure 5 are a result 
of the fact that the position at which the equilibrium concentration equals 
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the Initial concentration moves with time. This In turn results from the fact 
that the ends of the equilibrium profile are not perfectly symmetric with 
respect to the middle position of the column as can be seen In figure 5. 
However, when the gravitational equilibrium effect (product of g and the 
buoyancy factor) Is small, the profile Is more symmetric, and the Irregulari- 
ties disappear. This Is Illustrated by the points In figure 5 which represent 
G = 1 and a buoyancy factor of 100. However, the essential features In figure 
5 are Independent of G and the buoyancy factor; thus, the fraction of equilib- 
rium achieved Is essentially Independent of these two parameters: g and 

®1 • 


In any experimental Investigation, measurement of only the upper and lower 
concentrations Is all that would be necessary to establish the "equilibrium" 
values for the gravitational effect. Therefore, one can use the fraction of 
equilibrium achieved at the top and bottom of a column as an overall measure of 
the fraction of equilibrium achieved. To obtain a single value for the top and 
bottom, the product of G and the buoyancy factor can be assumed to be small. 
With these two assumptions, figure 7 Is plotted: the fraction of total change 

of concentration at the top (and bottom) Is the abscissa, and the log of time 
Is the ordinate. The dotted line represents the default values of the diffu- 
sion coefficient (10~ 5 cm 2 /sec) and height (1 cm). The time required to 
achieve over 96 percent equilibrium Is not graphed In figure 7 because It has 
nearly a vertical slope. Remember that equilibrium Is never achieved In a 
finite time frame. 


Effect of D and h 

Via the computer program In table III, the concentration function profiles 
(similar to those In figure 5) can be calculated for various values of D and 
h. The values of the concentration function at the ends (top and bottom) are 
again plotted as a function of time In figure 7. Curves for a constant value 
of D/h 2 are superlmposable. Undoubtedly, thlsfact could be proven mathe- 
matically starting with equation (28). 

From figure 7, the time, t, required to achieve a certain degree of equi- 
librium Is seen to be Inversely proportional to the square of h. Because of 
this, the curves In figure 7 can be combined Into the single curve In figure 8, 
where tO/h 2 Is plotted as a function of the fraction of equilibrium (concen- 
tration change) achieved. Via this graph and the values of D and h, one 
can readily determine the fraction of equilibrium achieved. If one Is content 
with 96 percent equilibrium, the following equation may be written as an 
approximation. 


t = 0.3 (33) 

where 0.3 Is a unitless value. For our default values, the time predicted 
by equation (33) Is 30 000 sec (about 9 hr). In contrast for sodium chloride 
In the ocean, the time required for 96 percent equilibrium Is predicted to be 
about 2xl0 16 sec (roughly a billion years). Because the time to approach 
equilibrium Is dependent only on D and h, and the D for all Ions are 
quite similar, roughly the same time to approach equilibrium In the ocean Is 
required by all salts. 
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Solid Solutions 


The default parameters for the computer program In table III represent a 
solute In a liquid solution because the diffusion coefficient used Is 10 - 5 
cm 2 /sec. The diffusion coefficient of the solute In a solid solution can be 
several orders of magnitude lower. Equation (33) would predict t to be 
900 hr when the diffusion coefficient of the solute Is 10 - 7 cm 2 /sec. This 
time Is likely too long for an ultracentrifuge experiment. 


Gases 

Gases exhibit diffusion coefficients In the range of 10 -4 cm 2 /sec, larger 
than those for solutes In liquid solutions. Thus according to equation (33), 
a 1 cm high column of gas would approach equilibrium In less than an hour. 

For a gas, the buoyancy factor Is Its molecular weight which Is In general 
larger than the buoyancy factor for a solute In a solution (Mi - \7ip). Thus, 
the gravitation - concentration effect should be more easily observable In a 
gas than In a liquid solution with the assumption that both are under the same 
gravitational conditions, e.g., 100 000 G with h = 1 cm. Under a lower 
gravitational force, h would have to be larger to observe the same difference 
In concentration between the top and bottom of the column. This would require 
a longer equilibrium time. If one sets the maximum time at 10 hr then equa- 
tion (33) gives a height of about 4 cm for the column of a gas. 


Summary of the Kinetics 

Figure 8 and equation (33) summarize our findings concerning the kinetics 
of gravitational Induced concentration gradients. The fraction of gravita- 
tional equilibrium achieved Is directly proportional to the diffusion coeffi- 
cient and Inversely proportional to the square of the height of the column. 

The time required to attain essentially 96 percent equilibrium Is given by 
equation (33). This equation shows that the diffusion Induced by gravity Is a 
very slow process; In general, very long times are required to approach equi- 
librium. Only for low values of h, about 1 cm or less. Is equilibrium 
approached In a reasonable time period, less than 10 hr. Gravitational - 
concentration equilibrium In the ocean would require a billion years. Gravi- 
tational equilibrium Is therefore not approachable because the convection In 
the ocean operates on a smaller time scale. 

Even If the kinetics are favorable, the equilibrium values themselves must 
be sufficiently different from the original concentration; see equation (10). 
This means that the gravitation constant must be large, of the magnitude 
achievable In an ultracentrifuge. 

Although attainment of concentration equilibrium Is a slow process, 
figure 5 Indicates that concentration equilibrium at the top and bottom of any 
column Is achieved Immediately. Unfortunately, this Is a very small gradient 
under normal gravity conditions: 0.001 percent/cm for sodium chloride In the 

ocean (or any column of water at 1 G). However, If a device were built that 
could measure such a small gradient, detection of the gravitational - concen- 
tration effect would cease to be a kinetics problem. 
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SPECIAL TOPICS 


Ionic Species 

As we have seen earlier, equation (14) describes the variation of the 
concentration of an Ionized salt In a solution subject to a gravitational 
field. However, equation (14) Is valid only when a single salt Is present In 
solution. One way to obtain a more general equation Is to start with the 
following equation. 

2( b ' a 1 ) - 0 (3 

where b^ Is the valence of an Ion (It Is negative for an anion); and the 
summation Is over all the cations and anions present. The equation merely 
states that an electrical charge gradient cannot exist, even In a gravita- 
tional field. Also required Is an equation that relates concentration gradi- 
ents to buoyancy factors. 


dC 1 

dh~ 

b 1 C 1 


dC 

dh 

Vj 



(35) 


This equation can be derived In the same manner as equation (7) by combining 
the equations for two different Ions so that the effect of electrical charge 
Is eliminated from the calculation. In equation (35) both subscripts, 1 and 
j, refer to any Ion present; note that the b's are negative for anions. By 
combining equation (34) and (35) followed by simplification, one obtains: 



(36) 


where the summations are again over all the Ions present; and j = - Bj 

= B s (1/b^ - 1/bj). When 1 and j are both anions or both cations, one may 
substitute Bi t |< - Bj |< for B^j where the subscript k represents any 
cation when 1* and j are both'anlons, and k represents any anion when 1 
and j are both cations. Thus, the use of equation (36) does not require 
knowledge of the values of B for Individual Ions, only those for salts. The 
Integrated form of equation (36) Is only an approximation because all the con- 
centrations, Cj, are functions of h. However, equation (36) allows one to 
demonstrate the general dependence of concentration ratios on the composition 
of a solution. 
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Our first example Is an aqueous solution of sodium chloride and ferric 
chloride. Figure 9 Is a plot of the concentration ratio of each Ion as a 
function of the relative concentration of the two salts: the cation fraction 

of the ferric Ion. The conditions assumed are 100 000 G's and h = 1 cm. The 
actual concentrations of the salts need not be specified. However, the plot 
In figure 9 Is more precise for very dilute solutions. The data used In the 
plot are available In table I. Note that the ratio NVN' for both cations 
varies with composition even though each of these cations can exist only as 
the chloride. Only when one salt predominate In a solution will equation (14) 
yield a precise ratio of NVN' for that salt. This Is Illustrated In 
figure 9 by the fact that the N"/N' of both the sodium and chloride Ions 
approach the value for sodium chloride (table I) near zero ferric Ion 
concentration. 

A more complex example Involves the mixture of two salts without a common 
Ion. This means that two cation and two anions are present, and thus the 
solution really contains four salts. The desired calculations can only be 
made when the buoyancy values of all four salts are known. The example to be 
used Is a solution of sodium chloride and potassium nitrate, which Includes 
following Ions: Na , K + , Cl” and NO3. The results of the calculations are 

plotted In figure 10. The same parameters are used as are used for figure 9: 
100 000 G's and h = 1 cm. Again one finds that the N"/N' ratio for each of 
the Ions varies with the concentration parameter. 

To plot the values of N"/N' as a function of all the Ions present In 
two salts (without a common Ion) would require at least a pseudo-three dimen- 
sional graph. Treatment of solutions with more than two salts would require 
multidimensional graphs. Thus, application of equation (36) to multlsalt sys- 
tems Is best accomplished by restricting calculations to the composltlon(s) of 
particular Interest. 


Osmosis 

Generally the van't Hoff equation (ref. 15) Is used to express osmotic 
pressure: P os = N2RT/V-) where 2 refers to the solute; and 1 to the sol- 

vent. Let us envision a vertical semlpermeable membrane separating a solution 
containing a solute from the pure solvent. A pressure equal to the osmotic 
pressure Is applied to the solution side so that the surfaces of both the pure 
solvent and the solution will be at the same level. Linder these conditions, 
there will be a hydrostatic pressure difference, AP^y, between the two sides 
of the membrane which Increases with depth according to the following equation. 


d AP 
dh 


hi _ 


= g(p - pj) 


( 37 ) 


where j refers to the pure solvent. For equilibrium to exist at all depths, 
the osmotic pressure would have to vary In the same way with h. To substan- 
tiate this, one can start with a different formulation of the osmotic pressure 
equation (ref. 15). 


(In N ) (in N ) 

P os = " RT ^7— = - RT “VT" 


’1 




(38) 
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where the final form of the equation uses our subscript notation. If Vj 
Is assumed to be Independent of h, differentiation of equation (38) yields: 


^os -RT d tn < N 1> 
dh = Vj dh 


(39) 


Substituting the value of d 4n(Nj)/dh from equation (9) Into equation (39) 
yields : 


dP 


os 

dh 


- V> 


V 


j 


g(p - pj) 


(40) 


The final form of the equation results from the fact that p\ = Mj/Vj. Com- 
parison of equations (37) and (40) shows the correspondence between the change 
In osmotic pressure with h and the change In hydrostatic pressure with h. 


Vapor Pressure 

Let us now apply equation (10) to the concept of vapor pressure. Assume 
that we have a pure liquid with a vapor that behaves Ideally. A liquid with a 
low vapor pressure will suffice. Actually, one could consider a solid just as 
well as a liquid. For such a material the fugacltles In equation (15) can be 
replaced by the vapor pressures because the fugaclty Is the Idealized vapor 
pressure. 

RT 8.n (jjr) = P( appl 1 ed ) V (41) 

where V Is the volume of a mole In the condensed state. Equation (41) states 
that the vapor pressure of a material Increases with an Increase In the applied 
pressure. However, generally the effect Is quite small. The effect of the 
earth's atmosphere on the vapor pressure of water at 20 °C Is only about 
0.07 percent. 

In equation (41) the applied pressure may be the hydrostatic pressure of 
the liquid Itself which changes with depth; P = pgh. Interpretation of equa- 
tion (41) In this manner tells us that the vapor pressure at the bottom of an 
Isothermal liquid In a gravitational field Is somewhat greater than that at 
the top. 

Equation (15) can be modified still further; the applied pressure can be 
assumed to be the vapor pressure. Thus, we have 

RT an (jt) - P"V (42) 

In this equation, P' Is the vapor pressure measured under vacuum (Langmuir) 
conditions while P" Is the value obtained under equilibrium (Knudsen) condi- 
tions. Equation (42) show us that from a theoretical viewpoint the two are not 
Identical. However, from a practical viewpoint, the difference Is negligible. 
For water which has a vapor pressure of about 17.5 mm of Hg at 20 °C, the 
Knudsen pressure would be only 0.0017 percent greater than the Langmuir value. 


18 



For ethanol which has a vapor pressure of about 42 mm of Hg at 20 degrees, the 
Knudsen pressure would be 0.008 percent greater than the Langmuir pressure. 

We have applied equations (41) and (42) to only pure materials. However, 
they may be applied just as validly to components In an Ideal solution. 


Two Continuous Phases In Contact 


If two phases are set up as side-by-side columns, and each phase Is 
assumed to be In equilibrium with respect to h, equation (18) holds for all 
the components of both phases. As a result 


f 

f 


II 


u 

I 

1,1 



(43) 


where the 1 and the 2 refer to the two different phases. The system Is 
most easily visualized In terms of a solid and a gas, but the concept Is valid 
for solid - solid, solid - liquid, liquid - liquid, and liquid - gas systems. 
Systems Involving two fluids requires the visualization of semlpermeable mem- 
brane between the phases. 


Let equilibrium be assumed to exist between the two phases at the higher 
altitude. In other words, f} = f 1 , 2 for each component. As a consequence 
of equation (43), f? i = fi,2* f° r each component. Thus, If two phases are 
In equilibrium with respect’to h, and In equilibrium with each other at one 
altitude, they are In equilibrium at all altitudes. Under these conditions no 
net transport of any species Is expected. Let us restate the conclusion for a 
more specific case. A column of a single phase solid setting In an atmosphere 
of Its own vapor Is not expected to be transported from a higher elevation to 
a lower elevation. This also means that extraction coefficients based on 
fugacltles do not vary with position In a gravitational field; this does not 
mean that an extraction coefficient based on concentrations will not vary 
with h. 


Two Phases In Contact, One Not Continuous 

Phase 2 Is now assumed to be noncontlnuous while phase 1 Is still con- 
tinuous. One portion of the noncontlnuous phase 2 Is allowed to contact the 
column of phase 1 at one altitude, and another portion of phase 2 Is allowed 
to contact the column of phase 1 at another altitude, a distance h below the 

first. In this case, equation (18) applies to the components In the continuous 

phase but not to those In phase 2. For phase 2 the fugaclty of a component Is 
determined by the pressure exerted by the continuous phase 1. 

RT 4n = ^ i,2f> = ^’ 2pi9h (44) 

where 2 the partial molal volume of species 1 In phase 2, and the 

pressure at the higher elevation Is assumed to be nil. 
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Next let us assume that equilibrium exists at the higher altitude. 


Combining equation (18) for phase 1 with equations (44) and (45) yields 

RT “(qi) ' <M ’ ' Vi , ’i )9h 

thus the fugacltles of species 1 are not the same In both phases at depth 
h. Transport of species 1 between the phases Is expected. If the buoyancy 
factor Is positive, transport at depth, h, will be from phase 1 to phase 2 
because the fugaclty Is greater In phase 1. This will cause a displacement of 
phase 1 from gravitational equilibrium; and eventually at the higher elevation, 
species 1 will tend to be transported from phase 2 to phase 1. Recall that 
such transport Is not predicted when both phases are continuous. The transport 
In the case of a noncontlnuous phase Is due to the fact that the lower portion 
of the noncontlnuous phase Is not subjected to a hydrostatic pressure equal to 
that of a column of phase 2, h In height. 

Under the conditions where transport does occur between different eleva- 
tions, the mechanism varies depending upon whether the phase 1 Is condensed or 
a gas. When phase 1 Is a gas, transport occurs by vaporization - condensation, 
and the rate of transport depends upon the vapor pressure of the species In 
phase 2 (or the vapor pressure of phase 2 Itself If It Is a pure material). 

Thus for very low vapor pressures, the transport could be undetectable. When 
phase 1 Is a liquid or solid, transport would occur by solution - diffusion - 
dissolution. The rate of transport will be highly dependent on the solubility 
of phase 2 In phase 1. Low solubility would mean that transport would essen- 
tially be nil. 

Transport can occur even If none of the material to be transported Is at 
the lower elevation; only nucleatlon sites are required. Consider the follow- 
ing example. A container Is filled with an aqueous saturated salt solution 
and a crystal of the salt Is adhering to the side and near the top of the con- 
tainer. Equation (46) predicts that the crystal of salt will tend to dissolve 
and recrystal 1 1 ze near the bottom of the container. (Evaporation of the water 
Is assumed to be nil .) 

When phase 2 Is a metal, transport can occur electrochemlcally If certain 
conditions are met. The metal at the upper position must be connected elec- 
trically to the metal In the lower position. The column of liquid must contain 
an electrolyte. Finally, the electrochemical transition between the metal and 
Its Ion must be reversible. Without the short circuit between the two Identi- 
cal pieces of the metal (each at a different elevation), a small voltage would 
exist between them. Such a voltage would not only be dependent upon the chem- 
ical potential difference between the two electrode but upon the transference 
numbers of the Ions Involved. Thus, a discussion of the voltage generated 
becomes too far outside the scope of this report. 


( 45 ) 


(46) 
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Multiphase Systems 


From the preceding sections, one can see that a multiphase system can be 
In equilibrium only when all the phases are continuous. If one phase Is non- 
contlnuous, complete equilibrium between various parts of such a phase Is not 
possible. This follows from equation (46). If diffusion between phases Is 
permitted, the system would be unstable as the example of a salt In a satu- 
rated water solution. In contrast If diffusion between phases Is not permit- 
ted, stability (but not equilibrium) of the system Is possible. In such a 
multiphase system, equilibrium with respect to h for continuous phases can 
still exist; equations (10) and (18) would apply to the components In these 
phases. Equilibrium with respect to h for noncontlnuous phases Is not pos- 
sible because N^| = for all the components of such phases. The fugacl- 
tles are determined by the hydrostatic pressure to which they are subjected; 
see equation (15). The hydrostatic pressure Is given by pgh, where the den- 
sity, p, Is that of a continuous phase surrounding the noncontlnuous phase. 


Solubility 

Consider again, two columns of different phases In a gravitational field, 
each phase being In equilibrium with respect to h. Let phase 2 be of a pure 
condensed phase, and phase 1 be a saturated solution of phase 2 In some sol- 
vent. Equation (10) Is applicable, and the N^'s can be replaced by S^'s, 
the solubilities. 


RT = 9hB2 ’ 1 (47) 

where B 2 j Is the buoyancy factor of phase 2 In phase 1; and S2,l Is the 
solubility of phase 2 In phase 1. 

Like equation (10) from which It was derived, equation (47) applies pre- 
cisely only to dilute solutions. Therefore, equation (47) applies best to 
relatively "Insoluble" materials. However, the partial molal volumes of such 
materials are sometimes difficult to measure. Therefore, the greatest appli- 
cation of equation (47) might be for more soluble materials where the equation 
would only be an approximation. For Instance, equation (47) predicts about a 
6 percent change (Increase In the downward direction) In the solubility of 
sodium chloride In water for h = 1 cm and 100 000 G's. The values of ^2,1 
and p used for the calculation are those In the vicinity of the saturation 
concentration. 

The ratio of solubilities of a noncontlnuous phase 2 In a continuous 
phase 1 may also be expressed as a function of h. The fugaclty ratio for 
phase 2 Is determined by the pressure exerted by phase 1. Thus, 

RT = V 2 P = V 2 pgh (48) 

where V 2 Is the molar volume of pure phase 2, and p Is the density of 
the saturated solution (phase 1). Because saturation Is assumed at both 
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altitudes, the fugaclty of the dissolved phase 2 Is equal to the fugaclty of 
solid phase 2. Thus, 


RT in 



v 2 pgh 


(49) 


Because Henry's law Is assumed to hold, the ratio f/N for an nonlonlzed mate- 
rial Is constant for any particular value of h. This constant may be evalu- 
ated from the equilibrium values of f and N. Therefore, the dependence of 
f/N on h Is the same as for the dependence of their equilibrium values on 
h, which can be ascertained by combining equations (10) and (18). 



where the i Is the partial molal volume of phase 2 In phase 1, not the 
molar volume of pure phase 2. By combining equations (49) and (50) and sub- 
stituting S's for N's because the concentrations correspond to solubilities, 

one obtains: 


« t - <>m v 2 - v 2>1 > <5,) 

The same equation may be derived for Ionized salts. Moreover, It may be 
derived In a manner which shows It Is not limited to dilute solutions. To 
apply equation (51) we again use a 1 cm high column of an aqueous solution of 
sodium chloride at 100 000 G's. The change In solubility predicted Is less 
than 2 percent, less than the change In solubility (6 percent) predicted when 
sodium chloride Is a continuous phase. 

The pgh In equation (51) may be replaced by P, the hydrostatic 
pressure. 


RT tn 




(52) 


From this equation one sees that the solubility change with h for a non- 
contlnuous phase can be attributed solely to the hydrostatic pressure of the 
solution. 


Equilibrium Constants 

The fact that solubility changes with h, leads to the question whether 
other types of equilibrium constants change with h. An equilibrium constant 
at the higher position In a gravitational field Is defined as follows. 
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K' 


(53) 



where 1 refers to the reactant species; j, to the product species; and the 
y's, to the coefficients of the species In the chemical reaction under consid- 
eration. By restricting the equilibrium constant to molecular (nonionic) 
reactions, equation (18) for each molecular species may be substituted Into 
equation (53). 

n ( f i >J ) ) 

"h y, )l 

But X(yiMi) = l( yjMj) so that 


K' 


The right side of equation (55) Is the definition of the equilibrium constant 
at the lower positions, K"; this shows that equilibrium constants Involving 
molecular species and based on fugacltles are Independent of depth In a gravity 
field. 


xp (K y^) - X(yjMj)} 


(54) 


3 ) 

n(q Vl ) 


(55) 



Polymers 

Let us write equation (10) for an addition polymer, assuming that \7 p 
= p' Vi where 1 represents the monomer; p, the polymer; and p ■ , the number 
of monomer units In the polymer. 


RT »n0 = p'(M 1 - V lP )gh (56) 

Combining this equation with equation (12) for the monomer gives the following 
relationship. 


N' “ N" 

P P 

Thus at gravitational equilibrium, the value of N p /N^P‘ Is a constant Inde- 
pendent of h, even though the polymer concentration may not be In equilibrium 
with the monomer concentration. Combining equation (57) with equation (12), 
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( 58 ) 


N" /N" \ 

jj? = ( ) ex Pt<P' " D9h("i " V>> 

Because Is expected to be greater than V^p, the ratio of N p /N^ should 

Increase with depth If equilibrium with respect to h exists. 

When equilibrium also exists between the polymer and the monomer concen- 
trations, equation (56) states that the equilibrium constant (based on con- 
centrations) Is Independent of h, providing Vn = P'^1- 1° a similar manner 

It can be shown that the equilibrium constant (in terms of concentrations) for 
any reaction Is Independent of h If the sum of the molal volumes of the 
products Is equal to the sum of the molal volumes of the reactants. 


Near the Critical Solution State 

The variation In concentration In a gravitational field seems to be most 
easily observed near a critical solution temperature. Most If not all Inves- 
tigations of this phenomenon have Involved two component systems. At the 
critical temperature and critical composition, such a system separates Into 
two phases. In the single phase region near the critical point, one might 
expect that molecules of like species would cluster together to form what 
might be called physical polymers (molecules held together by physical not 
chemical forces). Acceptance of this concept of physical polymers near the 
critical point would mean that equation (56) could apply. Via experimental 
data, equation (56) could be used to calculated a value for p ‘ ; p' would be 
an estimate of the number of molecules In a cluster (or polymer) of the 
component. 


Extrapolation to Higher Concentration 


To apply equation (7) to solutes of higher concentration, one must express 
the density In terms of the mole fraction of species 1. In this treatment the 
solution will be considered to contain only one solute. 


M 1 + N i (H 1 - Mj ) 
Vj * N,(V, - V 


(59) 


Substitution of this equation Into equation (7) and expansion of reciprocals 
to only the first power of N^ results In the following: 


1 

N 


1 



1 


dN 1 ‘ (M 1 " V 1 p j )g RT 


(60) 



where 1 refers to the solute; j, to the solvent; and Is the density 
of the solvent. To Integrate equation (60) one must still assume a constant 
and noncompresslblllty of the solution. 



where the subscript x refers to the value of n][/n{ calculated from 
equation (10). Equation (12) Is then used to calculate Nj. Finally, 
equation (61) can be used to obtain a closer approximation for the ratio 

n\/nJ . 


Suspensions 

Let us now return to Perrin's work (ref. 5) Involving a suspension of 
particles In a liquid. The fact that equation (5) Is obeyed leads to an 
Interesting Inference. Note that equation (5) treats each suspended particle 
as a molecule, and that equation (5) specifies a mole of suspended particles. 
This Is clear from the occurrence of Avogadro's number In the equation. If 
one attributes a molecular weight, M s , and a molecular volume, V s , to the 
particles, equation (5) becomes 

RT * n (!r) = (M S - V )9h < 62 > 

Because equation (62) has the same form as equation (10), and equation (10) 
can be viewed as a balance of chemical and gravitational potentials, one can 
postulate that a similar balance gives rise to equations (5) and (62). Thus, 
a chemical potential could be attributed to suspended particles In a fluid, 
just as one Is attributed to a solute In a solvent. However, remember that 
for the suspension, an Avogadro's number of particles Is taken as a mole of 
the material . 


Suggested Microgravity Experiment 

Equation (52) tells us that via pressure one can Increase the solubility 
of a solute. Releasing the pressure would cause crystals of the solute to 
form. The purity, size, and perfection of these crystals would depend on the 
slowness at which the pressure Is released. In 1 G, a difference In density 
between the crystals and the solution can Induce convection that can limit the 
perfection of the crystals. However, under microgravity conditions, the con- 
vection would be minimal, and the resultant crystals should be more perfect. 

An example would be the recrystalllzatlon of sodium chloride from water. 
The difference between the molar volume of the solid and Its molal volume In 
solution Is about 9.3 cc. With an applied pressure of 100 atm, an Increase In 
solubility of about 4 percent Is expected. For a liter of solution and a 
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solubility of 0.357 g/cc under zero applied .pressure, the expected yield Is 
14 g. This Is more than required to test the method. 


SUMMARY OF FINDINGS 


Under equilibrium conditions, the fugaclty of molecular and Ionic species 
Is a function of the position In a gravitational field: equation (18). This 

effect Is partly due to the hydrostatic pressure of the material above It: 
equation (15). The rest of the effect Is associated with a concentration gra- 
dient necessary for equilibrium: equation (10). For multiple Ionic species, 

equation (36) can be applied. The concentration effect Is generally less than 

0.1 percent for a meter column at 1 G. But under ultracentrifuge conditions 
of 100 000 G's, over a 10 percent difference Is generally expected between the 
concentration at the top of a 1 cm column as compared with Its concentration 
at the bottom, tables I and II. The ocean's depth could lead to even greater 
concentration gradients for the salts present. However, the existence of fluid 
convection Inhibits the attainment of such equilibrium. 


Modifications of the usual diffusion equations are necessary to treat the 
achievement of concentration equilibrium In a gravitational field: equations 

(27) and (28). The Important parameters for approaching equilibrium are the 
diffusion coefficient and the height of the column. Figures 7 and 8 give the 
time required to achieve various fractions of equilibrium. Equation (33) gives 
the time required to achieve 96 percent equilibrium. For 96 percent equilib- 
rium In a 1 cm high column, about 9 hr are required for a liquid solution. As 
the height Increases, the time required becomes extremely long; gravitational 
equilibrium In the ocean requires about a billion years. 


bulk mass trans- 
1s the existence 


Under certain conditions a gravity field can give rise to 
port: equation (46). Another effect of a gravitational field 

of solubility gradients: equations (47) and (51). However, equilibrium 

constants and extraction coefficients based on fugaclty are not affected b 
gravity field. Moreover, the dependence of fugaclty on depth and pressure 
helps explain some details of osmotic pressure and vapor pressure. 
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TABLE I. - PARTIAL MOLAL VOLUMES AND CONCENTRATION RATIOS OF SOLUTES IN AQUEOUS SOLUTION IN A GRAVITATIONAL FIELD 

(Temperature = 20 °C) 


Solute 

Species 

per 

molecule 

Mole 

weight 

Part, 
molal 
vol . 

Buoyancy 

factor 

Cone. 

ratio 

(a) 

Cone. 

ratio 

(b) 

Cone. 

ratio 

(0 

Cone. 

range 

percent 

Cone. 

ratio 

(d) 

Fuga. 

ratio 

(e) 

Acetic acid 

1 

60.05 

51.05 

8.955 

1.000 

1.036 

1.433 

30 

1.028 

1.273 

Acetone 

1 

58.05 

66.18 

-8.10 

.9999 

.9678 

.7215 

9.5 

.9706 

1.263 

Ammonium chloride 

2 

53.5 

36.40 

8.531 

1 

1.034 

1.409 

.5 

1.033 

1.240 

Ammonium hydroxide 

1 

35.05 

42.42 

-7.33 

.9999 

.9709 

.7444 

58 

.9799 

1.151 

Ammonium sulfate 

3 

132.1 

53.09 

26.32 

1 

1.111 

2.885 

.5 

1.109 

1.702 

Barium chloride 

3 

208.2 

25.79 

60.80 

1 

1.277 

11.55 

.5 

1.275 

2.312 

Calcium chloride 

3 

110.9 

20.34 

30.19 

1 

1.129 

3.371 

.5 

1.130 

1.563 

Cesium chloride 

2 

168.3 

41.37 

63.41 

1 

1.290 

12.84 

2 

1.290 

1.969 

Citric acid 

1 

192.1 

115.4 

76.57 

1 

1.361 

21.80 

1 

1.364 

2.166 

Cobalt chloride 

3 

129.8 

13.50 

38.77 

1 

1.168 

4.761 

.5 

1.168 

1.686 

Copper sulfate 

2 

159.6 

-2.36 

80.99 

1 

1.385 

26.04 

.5 

1.387 

1.901 

Creatinine 

1 

113.1 

86 

27.06 

1 

1.115 

2.972 

.5 

1.119 

1.576 

Dextran 

1 

72 000 

43 257 

28 699 

1.122 

(f) 

(f) 

5.5 

(f) 

(f) 

Ethanol 

1 

46.07 

55.30 

-9.17 

.9999 

.9637 

.6911 

9 

.9713 

1.203 

Ethylene glycol 

1 

62.07 

54.62 

7.427 

1 

1.030 

1.348 

9.5 

1.029 

1.283 

Ferric chloride 

4 

162.2 

23.30 

34.71 

1 

1.149 

4.044 

.5 

1.149 

1.921 

Formic acid 

1 

46.03 

33.63 

12.34 

1 

1.050 

1.643 

7 

1.047 

1.203 

Fructose 

1 

180.1 

110.1 

69.77 

1 

1.324 

16.58 

52 

1.259 

2.064 

Glucose 

1 

180.1 

111.8 

68.22 

1 

1.316 

15.58 I 

4.5 

1.310 

2.064 

Glycerine 

1 

92.09 

70.04 

21.96 

1 

1.092 

2.420 

58 

1.069 

1.448 

Hydrochloric acid 

2 

36.47 

18.28 

9.083 

1 

1.037 

1.441 

1.5 

1.036 

1.158 

Inul in 

1 

5 200 

3 227 

1 969 

1.007 

2 767 

(f) 

5.5 

2 344 

(f) 

Lactose 

1 

342.3 

205.6 

136.4 

1 

1.731 

242.6 

7.5 

1.694 

3.965 

Lead nitrate 

2 

331.2 

-2 460 

1 423 

1.005 

307.2 

(f) 

.5 

269.2 

3.793 

Lithium chloride 

2 

42.4 

17.87 

12.24 

1 

1.050 

1.637 

.5 

1.051 

1.186 

Magnesium chloride 

3 

95.23 

17.46 

25.91 

1 

1.109 

2.837 

.5 

1.109 

1.467 

Magnesium sulfate 

2 

120.3 

-1.78 

61.08 

1 

1.278 

11.69 

.5 

1.277 

1.623 

Maltose 

1 

342.2 

205.6 

136.4 

1 

1.731 

242.6 

.5 

1.752 

3.965 

Manganese sulfate 

2 

151 

3.741 

73.62 

1 

1.344 

19.36 

.5 

1.342 

1.836 

Mannitol 

1 

182.1 

116.7 

65.35 

1 

1.300 

13.88 

.5 

1.290 

2.081 

Methanol 

1 

32.04 

37.81 

-5.74 

.9999 

.9771 

.7936 

56 

.9831 

1.137 

Nickel sulfate 

2 

154.7 

-11.4 

83.13 

1 

1.397 

28.39 

.5 

1.390 

1.864 

Nitric acid 

2 

63.02 

29.08 

16.94 

1 

1.070 

1.978 

1.5 

1.070 

1.288 

Phosphoric acid 

1 

98 

43.27 

54.66 

1 

1.246 

9.027 

.5 

1.240 

1.483 

Potassium bromide 

2 

119 

33.63 

42.65 

1 

1.187 

5.567 

1 

1.185 

1.614 




TABLE I. - Concluded 


Solute 

Species 

Mole 

Part. 

Buoyancy 

Cone. 

Cone. 

Cone. 

Cone. 

Cone. 

S91 


per 

weight 

mol al 

factor 

ratio 

ratio 

ratio 

range 

ratio 



molecule 

vol . 


(a) 

(b) 

(0 

percent 

(d) 

(e) 

Potassium carbonate 

3 

138.2 

14.37 

41.26 

1 

1.180 

5.263 

.5 

1.179 

1.744 

Potassium chloride 

2 

74.55 

26.99 

23.75 

1 


2.601 

1.5 

1.099 

1.349 

Potassium chromate 

3 

194.2 

39.45 

51.55 

1 


7.964 

.5 

1.229 

2.185 

Potassium dichromate 

3 

294.2 

88.98 

68.35 

1 

1.316 

15.66 

1 

1.313 

3.267 

Potassium iodide 

2 

166 

46.91 

59.51 

1 

Bmul 

10.97 

. 5 

1.273 

1.950 

Potassium nitrate 

2 

101.1 

38.61 

31.21 

1 j 

KSB9 

3.512 

1 

1.134 

1.502 

Potassium oxalate 

3 

166.2 

46.96 

39.71 

1 

1.173 

4.946 

1 

1.173 

1.952 

Potassium dihydrogen 

2 

136 

38.45 

48.78 

1 

1.216 

7.123 

1 

1.214 

1.729 

phosphate . 

Potassium monohydrogen 

3 

174.1 

25.02 

49.69 

1 

1.221 

7.391 

.5 

1.217 


phosphate 
Potassium sulfate 

3 

174.2 

35.40 

46.26 

1 

1.204 

6.436 

.5 

1.203 

2.016 

Propylene glycol 

1 

76.09 

70 

6.070 

1 


1.276 

9.5 

mmwm 

1 . 358 

Silver nitrate 

2 

169.8 

21.04 


1 

1.349 

19.97 

. 5 

1.J42 

1.981 

Sodium acetate 

2 

82.04 

39.48 

21.24 

1 

1.089 

2.352 

.5 

1.087 

1.391 

Sodium bicarbonate 

2 

84.02 

23.74 

30.11 

1 

1.128 

3.360 

1 


1.402 

Sodium bromide 

2 

102.9 

22.95 

39.95 

1 

1.174 

4.993 

.5 

1.173 

1.513 . 

Sodium carbonate 

3 

106 

-3.66 

36.55 

1 

1.158 

4.355 

.5 

1.157 

1.532 

Sodium chloride 

2 

58.45 

17.67 

20.37 

1 


2.270 

.5 

1.086 

1.265 

Sodium hydroxide 

2 

40.01 

-4.55 


1 

1.093 

2.452 

.5 


1.174 

Sodium molybdate 

3 

205.9 

29.58 

58.76 

1 

1.266 

10.64 

.5 

1.265 

2.290 

Sodium nitrate 

. 2 

85.01 

27.39 


1 

1.122 

3.185 

.5 

1.120 

1.407 

Sodium monohydrogen 

3 

141.9 

.7063 

■ 47.09 

1 

1.208 

6.654 

.5 

1.205 

1.770 

phosphate 
Sodium sulfate 

3 

142 

14.77 

42.41 

1 

1.186 

5.513 

.5 

1.187 

1.771 

Sodium tartrate 

3 

194 

58.69 

45.09 

1 

1.198 

6.140 

1 

1.199 

2.183 

Sodium thiosulfate 

3 

158.1 

28.99 


1 

1.189 

5.650 

.5 

1.190 

1.889 

Sodium tungstate 

3 

293.9 

30.57 

87.75 

1 

1.423 

34.19 

.5 

1.426 

3.264 

Strontium chloride 

3 

158.5 

19.63 

46.28 

1 

1.204 

6.443 

. 5 

1.206 

1.892 

Sucrose 

1 

342.3 

209.3 

132.5 

1 

1.705 

207.6 

58 

1.530 

3.965 

Sulfuric acid 

2 

98.08 

31.83 

33.06 

1 

1.142 

3.784 

1 

1.139 

1.484 

Trichloroacetic acid 

2 

163.4 

81.90 

40.69 

1 

1.177 

5.144 

1 

1.179 

1.930 

Urea 

1 

60.06 

43.26 

16.76 

1 

1.069 

1.963 

1 

1.067 

1.273 

Zinc sulfate 

2 

161.4 

-5.58 

83.52 

1 

1.399 

28.83 

.5 

1.397 

1.915 


a Dilute solution; 1 G; h = 1 m. 

^Dilute solution; 100 000 G's; h = 1 cm. 
c Dilute solution; 1 G; h = 10 000 m. 

^Conc. indicated in 'range 1 column; 100 000 G's; h = 1 cm. 
^100 000 G's; h = 1 cm. 

'Exceeds 10 to the 38th power. 





















TABLE II. - PARTIAL MOLAR VOLUMES AND CONCENTRATION RATIOS OF SOLUTES IN METALS AND ORGANIC SOLVENTS 


IN A GRAVITATIONAL FIELD (Temperature = 20 °C unless otherwise specified) 


Solute/solvent 

Mole 

weight 

Part, 
molal 
vol . 

Buoyancy 

factor 

Cone. 

ratio 

(a) 

Cone. 

ratio 

(b) 

Cone. 

ratio 

(c) 

Cone. 

range 

percent 

Cone. 

ratio 

(d) 

Fuga. 
rati o 
(e) 

Benzene 

78.11 

97.86 

-122 

0.9995 

0.6105 

0.0072 

10 

0.6595 

1.369 

/stannic chloride 










Bromocamphor 

231.1 

211.9 

-83.4 

.9996 

.7130 

.0339 

11.12 

.7250 

2.551 

/tristearin; 18 *C 










Ethanol/acetone 

46.07 

60.48 

-1.91 

.9999 

.9923 

.9258 

5.455 

.9957 

1.203 

Sodium nitrate 

84.99 

-1689 

1241 

1.004 

136 

(f) 

.07 

42.41 

1.399 

/methylamine; 25 °C 









Stannic chloride 

260.5 

161.2 

115.6 

1 

1.592 

105.1 

10 

1.664 

2.853 

/benzene 










Tristearin 

891.4 

511.2 

365.3 

1.001 

4.395 

(f) 

20 

3.883 

37.06 

/bromocamphor; 18 °C 









Tristearin 

891.4 

505.5 

364.4 

1.001 

4.827 

(f) 

20 

4.090 

47.03 

/bromocamphor; 0 °C 









Tri stearin 

891.4 

564.6 

398.7 

1.001 

3.652 

(f) 

20 

3.113 

18.10 

/bromocamphor; 90 °C 









Aluminum/copper 

63.54 

6.633 

44.36 

1 

1.195 

5.964 

18 

1.196 

1.291 

Aluminum/magnesium 

26.98 

15.50 

2.772 

1 

1 

1 

2 

1.015 

1.114 

Aluminum/zinc 

26.98 

10.34 

-41.6 

.9998 

.8457 

.1873 

9.59 

.8649 

1.114 

Copper/aluminum 

63.54 

6.634 

44.37 

1 

1.195 

5.965 

18 

1.196 

1.291 

Copper/gold; g 

63.54 

8.442 

-93.4 

.9996 

.6863 

.0232 

5 

.7028 

1 1.291 

Iron/ aluminum 

55.84 

-.529 

57.28 

1 

1.259 

10.03 

1 

1.266 

1 1.252 

Iron/copper 

55.84 

6.260 

1.267 

1 

1 

1 

.96 

.9901 

1.252 

Lead/gold; g 

207.1 

20.01 

-171 

.9993 

.5020 

.0010 

5 

.4613 

2.302 

Lead/mercury; g 

207.1 

18.02 

-35.6 

.9998 

.8663 

.2380 

10 

.8586 

2.302 

Lithium/mercury; g 

6.939 

5.912 

-65 

.9997 

.7697 

.0730 

2.5 

.7465 

1.028 

Magnesium/aluminum 

24.31 

13.01 

-10.7 

.9999 

.9575 

.6483 

1 

.9785 

1.102 

Manganese/aluminum 

54.93 

8.328 

32.40 

1 

1.139 

3.684 

.5 

1.254 

1.247 

Nickel/aluminum 

58.71 

7.349 

38.78 

1 

1.168 

4.764 

1 

1.228 

1.266 

Oxygen/copper 

16 

5.816 

-35.8 

.9998 

.8655 

.2359 

.015 

2.747 

1.066 

Phosphorus/copper 

123.8 

-83.1 

866.7 

1.003 

32.73 

(f) 

.04 

4.941 

1.646 

Potassium/mercury; g 

39.10 

5.504 

-34 

.9998 

.8717 

.2535 

5 

.8818 

1.170 

Silicon/aluminum 

28.08 

13.95 

-9.10 

.9999 

.9640 

.6930 

8 

.9722 

1.119 

Silver/gold; g 

107.8 

15.69 

i -182 

.9992 

.4797 

.0006 

5 

.4849 

1.543 

Sodium/mercury; g 

22.98 

3.236 

1 -20 

.9999 

.9224 

.4462 

5 

.9287 

1.096 

Tin/gold; g 

118.6 

21.97 

-281 

.9988 

.3224 

IE— 5 

5 

.3569 

1.612 

Zinc/copper 

65.37 

21.90 

-118 

.9995 

.6217 

.0086 

7 

.6776 

1.300 


^Dilute solution; 1 G; h = 1 m. 
b Dilute solution; 100 000 G's; h = 1 cm. 

^Dilute solution; 1 000 000 G's; h = 1 cm. 

d Conc. indicated in 'range' column; 100 000 G's; h = 1 cm. 

^100 000 G's; h = 1 cm. 

'Exceeds 10 to the 38th power. 

^Obtained from graph. 












TABLE III. - KINETICS PROGRAM FOR DIFFUSION IN A GRAVITATIONAL FIELD 

10 REM GRAV/DIF 

11 GOTO 101: REM SET UP PARAMETERS 

28 REM LINES 28-31 CALC FLOW; AD=MAX 

29 T = T + U: FOR X = 2 TO N:S2 = F(X):S1 = C(X - 1,G) - C(X,G):F(X) = D * (SI/ 

L + SQR ((1 + C(X t G) ) * (1 + C(X - 1 , G) ) ) * Y) : IF T = U THEN 31 

30 IF ABS ( F(X) - S2) > AD(X) THEN AD(X) = ABS (F(X) - S2):BD(X) = T 

31 NEXT X 

32 FOR X = 1 TO N:S2 = C(X,G) : C ( X f G) = (F(X) - F(X + 1 ) ) / L * U + C(X,G): IF 
T = U THEN 35 

33 IF ABS ( C( X, G) - S2) > AC(X) THEN AC(X) = ABS (C(X,G) - S2):BC(X) = T 

34 S3 = ABS ( C( X , G) - C(X - 1 ,G) - S2 + S1):S1 = S2: IF S3 > AE (X) THEN 

AE ( X ) = S3 : BE(X) = T 

35 NEXT X: REM LINES 32-35 CALC NEW CONCENTRATIONS 

36 IF T < .9999999 * S * U THEN 29: REM SELECTS VALUES TO BE PRINTED AS 

DETERMINED BY DATA SET 

38 HOME : PRINT S$: PRINT "GRAD. T="; INT(T) : FOR X = 2 TO N:W = (C(X,G) - 

C(X - 1 , G) ) / ( B( X) - B ( X - 1)) * (1 + B( X) ) / (1 + C(X , G) ) * 1000: IF X / 

ST < > INT (X/ST) THEN 40 

39 IF W < > 0 THEN PRINT ■ ";X;")";.l * INT (W); 

40 WW = WW + W: NEXT X:WW * UW / (N - 1): PRINT " AVER) “ ; . 1 * INT (WW):WW = 0 

41 PRINT "CONC. T=" ; INT (T): FOR X = 1 TO N STEP ST:W = C(X,G) / B( X) * 

1000 

42 IF W < > 0 THEN PRINT 11 " ; X 1 * INT (W); 

43 NEXT X: PRINT : PRINT E$ 

44 G(G) = T:G = G + 1 

45 FOR X = 1 TO N:C(X,G) = C(X,G - 1): NEXT X 

59 READ S 

60 IF G < 6 OR ABS (C(N / 2 + 1,G)) < .95 * ABS (B(N / 2 4 - 1)) THEN 29 

70 PRINT S$: PRINT 

71 IF XX$ < > 11 Y" THEN 85: REM BYPASSES ACTUAL GRAD AND CONC DATA PRINT OUT 

72 FOR G1 = 1 TO G - 1 

73 HOME : PRINT "GRAD. T="; INT (G(G1)): FOR X = 2 TO N:W = (C(X,G1) - 

C(X - 1 , G1 ) ) / L: IF X / ST < > INT (X / ST) THEN 75 

74 IF W < > 0 THEN PRINT " 11 ; X ; " ) " ; W ; 

75 NEXT X: FOR X = 2 TO N:WW = WW + W: NEXT X: PRINT " AVER)";W / (N - 1):WW 
= 0 

77 PRINT "CONC. — T="; INT (G(G1 ) ) : FOR X = 1 TO N STEP ST:W = 1 + C(X,G1): 
PRINT " " ;X; " ) " ;W; : NEXT X: PRINT 

81 NEXT G1 

85 PRINT : PRINT "TIME OF MAX FLOW" 

90 FOR X = 2 TO N: PRINT " " ;X; ") 11 ;BD(X) ; : NEXT X: PRINT : PRINT "TIME OF MAX 
CONC CHANGE": FOR X = 1 TO N: PRINT " " ;X; 11 ) " ;BC(X) ; : NEXT X: PRINT : PRINT 
"TIME OF MAX CONC GRAD" 

91 FOR X = 2 TO N: PRINT " " ; X ; " ) " ; B E ( X ) ; : NEXT X: PRINT 

95 PRINT E$ 

100 HOME : PRINT "FINISHED": END 

101 HOME : REM EVERY QUESTION HAS A DEFAULT VALUE 

103 PRINT "ENTER DIFFUSION COEFFICIENT IN CM2/SEC" : INPUT D$:D = VAL (D$): IF 
D = 0 THEN D = IE - 5 



TABLE III. - Concluded. 


105 PRINT : PRINT "ENTER THE NUMBER OF G's": INPUT G2$:G2 = VAL (G2$): IF G2 
= 0 THEN G2 = 1E5 

107 PRINT : PRINT "ENTER BUOYANCY FACTOR IN GRAMS/MOLE": INPUT BU$:BU = VAL 
(BU$) : IF BU = 0 THEN BU = 100 

109 PRINT : PRINT "ENTER HEIGHT IN CM": INPUT H$:H = VAL (H$): IF H = 0 THEN 
H = 1 

111 PRINT : PRINT "ENTER TEMPERATURE IN CENTIGRADE" :H$ = " INPUT H$:TEMP = 


VAL (H$) : IF TEMP = 0 THEN TEMP = 293 

112 PRINT : PRINT "ENTER NUMBER OF INCREMENTS": INPUT N$:N = VAL (N$): IF N = 

113 PRINT : PRINT "DO YOU WANT A PRINT OUT? (Y/N)": INPUT H$: IF H$ < > "N“ 
THEN S$ = CHR$ (4) + "PR#1":E$ = CHR$ (4) + "PR#0" 

114 PRINT : PRINT : "DO YOU WANT ACTUAL CONC AND GRAD DATA PRINTED OUT? 

INPUT XX$ 

115 PRINT : PRINT "ENTER STEP INCREMENT IN PRINT OUT: 1, 2, 3, OR 4":INPUT 
ST$:ST = VAL (ST$): IF ST = 1 OR ST = 3 OR ST = 4 THEN 117 


116 = 2 

117 HOME : READ S:A = 1: DIM C(N,15): DIM F(N + 1): DIM B(N): DIM G( 15): DIM 
AD(N) : DIM BD(N) : DIM AC(N): DIM BC(N) 

118 G = 1: DIM AE(N) : DIM BE(N) 

119 Y = BU * G2 * 980 / 8.31E7 / TEMP: L =H/N:U=L*L/2/D/A 

120 U1 = INT ( LOG (U) / LOG (10) ) :U2 = 10 * Ul: IF U > 5 * U2 THEN U = 5 * 


U2 : GOTO 123 

121 IF U > 2 * U2 THEN U = 2 * U2: GOTO 123 

122 U = U2 

123 FOR X = 1 TO N:Q = Q + EXP (Y * ((X - .5) * L - .5 * H)): NEXT X 

124 0 = N / Q:Q = LOG (Q) 

125 FOR X = 1 TO N:B(X) = - 1 + EXP (Y * ((X - . 5) * L - . 5 * H) + Q) : IF 
B(X) = 0 THEN B ( X ) = - IE - 38 

127 NEXT X: REM LINES 120-126 CALCULATE EQUILIBRIUM CONCENTRATIONS, B 

129 PRINT S$: PRINT "BUOY=";BU;" HEIGHT=" ;H; " # OF SLICES- " ;N; " DIFF 
COEF=“:D:" ":G2;"G's TIME INCREMENT^" ;U: PRINT 

130 PRINT "EQUIL VAL FOR X = 1 TO N: PRINT " " ;X; " ) " ;B(X) ; : NEXT X: PRINT 
: PRINT E$ 

140 GOTO 29 „„„„ 

1000 DATA 1,2,5,10,20,50,100,200,500,1000,2000,5000,10000 
1002 DATA 20000,50000,100000,200000,500000 
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Figure 5. - Profiles of attainment of concentration equilib- 
rium in gravitational field (time on curves is in seconds; 
h = 1 cm ; D = 10' 5 cm 2 /sec). 
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Figure 6. - Profiles of attainment of concentration gradient 
equilibrium in gravitational field {time on curves is in 
seconds; h - 1 cm ; D - 10 -5 cm 2 /sec). 
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